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o : Abstract 

Confining a quantum particle in a compact subinterval of the real line with Dirich- 
£f) • let boundary conditions, we identify the connection of the one-dimensional fractional 

Schrodinger operator with the truncated Toeplitz matrices. We determine the asymptotic 
r*j | behaviour of the product of eigenvalues for the a-stable symmetric laws by employing 

the Szego's strong limit theorem. The results of the present work can be applied to a 
recently proposed model for a particle hopping on a bounded interval in one dimension 
whose hopping probability is given a discrete representation of the fractional Laplacian. 
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values. 
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O ■ 1 Introduction 

q | It is well known that a stable law [TJ El El HI [5] , which is a direct generalization of the Gaussian 

distribution, is generated by the parameters (a,(3,c,r) where: 
• i . 

— a is the characteristic exponent that determines the degree of leptokurtosis and the fatness 
& ■ of the tails. In the present work we consider a G (1,2). 

— /3 is the skewness parameter which characterizes the degree of asymmetry of the Levy mea- 
sure and takes values in the interval [—1,1]. 

— c is the scale parameter with range (0, oo) and measures scale in place of standard deviation. 

— t is the location parameter which saturates the set of real numbers and shifts the distribution 
to the left or right. 

It is often denoted by S a (j3, c, r). We adopt this notation to declare different parametric families 
of operators one might consider at quantum level. Our study will be focused on S a (c) (a-stable 
symmetric) operators. 

In the literature the papers P, [7j claim that the time-independent fractional Schrodinger 
equation with infinite one-dimensional square well potential has the same eigenf unctions as for 
the standard non- fractional case, only with modified energies. Recently in [8], it was argued via 
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a proof by contradiction, that the ground state cannot be a solution either in the interval a G 
( — 1, 1) (unless a = 0) or in (1, 2). Nevertheless, by using the Griinwald-Letnikov definition for 
the fractional derivative, one can numerically evaluate the eigenvalues and extract information 
about the asymptotic behaviour of the product of eigenvalues without knowing their explicit 
form. To prove this statement we organize our paper as follows. 

In Section II the square well potential with perfectly rigid walls serves as a simpified model 
from which one can derive and compare quantities (such as the determinant and trace) stem- 
ming from analytical and discretization methods. At discretized level one encounters the real 
symmetric Toeplitz matrix (TlOl) which forms a subset of the class of symmetric centrosym- 
metric matrices [13]. The eigenvalue problem can be solved exactly (fTTl) and the asymptotic 
behaviour of the product of eigenvalues is an easily accessible task since a closed expression for 
the determinant of this matrix can be determined. 

In Section III the situation for time-independent fractional operators in the class S a (/3,c), 
is more involved. Adopting the Griinwald-Letnikov definition for fractional derivatives with 
anisotropic coefficients, retaining Dirichlet boundary conditions and using the discretization 
method, we end up with the finite dimensional matrix ( 1241) . This matrix is recognized to be a 
Toeplitz matrix with well known properties. Restricting to the subclass of S a (c) laws one can 
prove that its eigenvalues are simple and non-negative thus positive definite. 

In Section IV the symbol (or generating function) / of the Toeplitz matrix (I24I) is shown 
to belong to the intersection of the Wiener and Besov spaces, and all requirements of Szego's 
strong limit theorem are satisfied except that the zeros of / are located at even integer multiples 
of 7r. Excluding these points from the complex unit circle the logarithm of the symbol can be 
expanded in Fourier modes while the determinant of fpHl) diverges as one might expect since 
the point spectrum of the associated operator is unbounded from above. Our final result for 
the asymptotic behaviour of the product of eigenvalues is captured by (jUJ) and fl5^|) which hold 
for every rational a in the superdiffusion region (1, 2). 

2 The square well potential as a toy model 

We study the asymptotic behaviour of the spectrum for a particle confined in an one-dimensional 
square well potential of length L with perfectly rigid walls. The potential is defined as 



The usual time-idependent Schrodinger equation with Dirichlet boundary conditions is 




(1) 



2m dx 2 



ip(x) 



Eip(x), < x < L 
0, Vx £ (0, L) 



(2) 



and has solutions the normalized eigenfunctions, 



The operator — 




(3) 




(4) 
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The product of N eigenvalues is 



(AT+1) 



2JV+1 



(5) 



where the Stirling's formula N\ w ^/2n(N + l^+ae ( N+1 ) for large iV has been applied. The 
sum of the N eigenvalues is 



N 



71=1 



/ h 2 TT 2 

\2mL 2 



N 



n 



71=1 



1 „ 1 



-iV 2 + -iV + - ) N. 
3 2 6 1 



(6) 



\2mL 2 

In the discretized method we consider a grid of N'— ordered points of the interval [0, L] 

x = < X\ < ■ ■ ■ < xn' < xn'+i = L (7) 

with equal spacing — Xi = e > 0, Vz = 0, • • • , N'. Then by employing the centered second 
difference estimator of the second ordered derivative, namely 



d 2 ,. , 4>(x + e) - 2ip(x) + 4>(x - e) 

— i>(x) = hm 

dx 2 e^o+ e 2 



(8) 



and denoting the values of the function ip at the grid points by ipi = ip(xi), the problem (j2J) is 
equivalent to 



(A N >) ij il>i = Eii()iJ = l,--- ,N' 
with i/;_ k =-^0 = = ipN'+k, V/c G N 



(9) 



where 



/ 2 -1 



.4 



h 2 



N' 



2me 2 



B N , 



h 2 



2me 2 



\ 




V 





The eigenvalues and eigenvectors of Bn> are [9] 

Afc = 2 | 1 + cos 



-1 2 



(10) 



/ N'xN' 



sm 



kir 
N' + l 
k(m + l)7r 
N' + l 



1 < k < N' 
< m < N' - 1. 



The trace of Aw is 
Tr(A N ') 



/i 2 



2me 2 

h 2 

2mL 2 



Tr(^ 



/, 2 



2mL 2 
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\ k=l V 



(N'+ 1) 2 2N 



(12) 
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while its determinant is given recursively by 



h 2 \ , / h 2 ^ 2 



det Ajv ' = 2 ( w j det An '~ 1 ~ { w ' det ^'~ 2 (13) 



with solution 

A 1 



The determinant could had equally been evaluated using the eigenvalues (TTTT) 

dets^ = n Afc = ri( 2 + 2c ° s ( ^ttt ) } = ^ + l (i5) 

fc=l fc=l ^ 



Increasing the number N' of intermediate points we approach the exact eigenvalues of problem 
with higher accuracy. 



3 The discretized fractional Schrodinger operator versus 
Toeplitz matrices 

We consider the fractional operator [10] 



1 / d a d a \ 

A(a,/3) = K+— + K- " , (16) 

V 7 cos(^) V ^ Q d(-x) a J v ; 

where i ; T ± are anisoptropic constants satisfying K > and K + + > 0. Operator (TT6|) can 
be casted into the equivalent form 

K f d a d a \ 

A(a,/3) = g-^r (1 + h (1 - /?) , " r (17) 

V ' 2cos(^) V dx a v HJ d{-x) a J K ' 

where [K a ] = ^ M |^ji — and 

K a = K + + K~, «G(0,2)/{1}, P= xl + K- > ( 18 ) 



For the case a = 1 one should consider an operator as the one proposed in [TO]. We study only 
the superdiffusive region 1 < a < 2 since the subdiffusive can be treated on equal footing with 
minor modifications. 

Definition 3.1. The derivatives in [Fj\ ) for LJ n functions are defined by 

f-1 

(^a \ (A a, 0) (x) 1 
- — -0 I (x) = lim — = lim — > w n {a)ib(x — (n — l)e), for x > 
7 n=0 

( d a \ (A a if) (x) 1 e 

( — r — -0 I (x) = lim - — — — = lim — > w n (ot)ib(x + (n — l)e), for x < L(19) 



and are called the right-handed (left-handed) Griinwald-Letnikov fractional derivatives [77] / 0. 
The weights w n (a) are given recursively by a 



w {a) 
w n (a) 



a + 1 



n 



w n -i(a) 



-l) n 



a 
n 



a (a — l)r(n — a) 
r(2-a)r(n + 1) 



Vn e N. 



(20) 



When a = [a], a positive integer, then iy n ([a]) = 0, Wn > [a] + 1. The only negative weight 
in the superdiffusive region is Wi(a) and the infinite series of weights converges to zero (see 
( lA.ip for the proof) 



n=0 



n: 



0. 



(21) 



Substituting a = 2 and /? = into f[P51) we recover the centered second difference estimator 
since w n {2) = 0, Wn > 3. 
The eigenvalue problem 



A(a, P)ip) (x) = \ip(x), ip e L 2 {0:L) 
with ip(x) = 0, Wxi (0,L) 
using a grid of n-ordered points, as in Section II, gives the discretized problem 



n 



with ip. 



^o = = ^ n+fc , Vfc G N. 



(22) 



(23) 



The approximation of the operator A(a,/3) by the sum of two matrices corresponding to the 
left- and right-fractional derivatives, results in a ~R nxn matrix with elements 



2 cos (*p) 



x < 



2wi 




for 


J 


= *; 




i = 


I,-- 


,n 


(1 + 


P)w Q +{l-p)w* 


for 


J 


— i - 


Hi; 


i = 


I,-" 


,n 


(1- 




for 


J 


= i - 


H2; 


i = 


1,... 


,n 


(1- 


(3)w A 


for 


J 


= i - 


H3; 


i = 


1,... 


,n 


(1- 


0)w n -! 


for 


3 


— i - 




i = 


1,2 




(1- 


p)w n 


for 


j 


— i - 


hn-1; 


i = 


1. 





1 

2 

3 



(24) 



Note that we do not allow the end points xq = and x n+ \ = L to be present in the calculation. 
The elements of the matrix parallel to the main diagonal are equal and (124j) is recognized to 
be a Toeplitz matrix. Renaming the entries of the matrix according to 

K„ 



(A(a,/3)L 



i.i 



2cos(^) e a 



n 



(25) 



Actually this is a variant of the Griinwald-Letnikov fractional derivative, in which the function evaluations 
are shifted to the right. 

2 Another way to compute the coefficients w n (a) is by their generating function w(z) = (1 — z) a which has 



the Taylor expansion w(z) = X)fc!lo(~ -0 



for \z\ < 1 and a £ [1,2]. 



In the subdiffusive region all weights are negative except wq- 
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its symbol in general is given by 



f(a,(3;z;n) = ^ fk(a,(3); 



z g 



(26) 



k=—n 



This is a rational function of z which is analytic everywhere on the complex plane except at 
the n-poles. On the counterclockwise complex unit circle T = {z G C : \z\ = 1} the generating 
function (126]) . for a large ordered sequence {x n }„ g N of intermediate points of the interval (0, L), 
is given by the expression 



lim f(a,/3;T;n) 

n— »oo 

u(a\ 9; oo) 



/(a, j3\ e ld ; oo) = u(a; 9; oo) + iv(a, {3; 9; oo), where 



(27) 



|cos(^)| e c 
K„ 



cos 



v(a, P; 9; oo) 



K a j3 
| COS (£2 )| e a 



a — cos # 
2 sin 



^ cos (^)___ 

fe=i 



a(a — 1) 



a) 



cos 



7T 



(28) 



. a(a - 1) ^ . T(fc+ 1 

sm9 — — r > sin w — — — — 

r(2-a^^ v 7 TYfca. 



cos 



\ 2 ) 



2 sin 



fe=i 



r(^ + 2) 



7T 



a 
2/ J 



(29) 



The infinite cosine series is non-negative definite for 9 ^ 2m7r, 27r-periodic and conveges to 
([28]) (see (IA.2D for the proof) for mT < 9 < 2n + mT, where T is its period. Function (1251) is 
periodic only for rational values of a since the sine and cosine functions have periods T\ = 2irrii, 
T 2 = ttt^- and therefore the product function is periodic only for 71% = ri\ (l — -). For example 
the value a = § gives (ni, ^2) = (4, 1) and T = 87r as it is also justified in the following plot. 









5 iV\ 


15 


/V 25 


30 35 



Figure 1: The plot of the function g{9) = |2 sin (§) \° cos [{9 — tt) (l — f )] for G [0, llvr]. 
We can shift the origin by substituting 9 — > 6* — 7r in (1251) and (129]) thus obtaining 



wfa;; 0: 00) 



f (a, 0; 00) 



|cos(^)| e c 
K a (3 



, 

I - 



, 
I - 



cos 



2/J 



sin 



0(1-- 
2 



I cos )| e c 

where — tt + mT < # < ir + mT. Note that for a = 2 and /3 = the symbol becomes 



/(2, /?; T; n) = u(2, T; n) = 2K a (l - cos 0) = K a 2 sin 







(30) 
(31) 

(32) 



with range [0, AK a ]. 

In the complex plane the symbol is graphically represented by a simple closed contour C. 




Figure 2: Symbol curves f(a, (3\ T; n) in the complex plane for the generating function ( 1271 ) 
with parameters a = 1.5, (3 = 0.8 and a = 1.2, (3 = 0.2 (the far right contour). The number of 
intermediate points is n = 200. 

For a given value of 9 e (0, 2%) the point (w, v) G C satisfies the equation 

-mX + (^mX = Asin^r. (33) 



The matrix fT24"j) contains the following subclasses of matrices 



{%) The matrix A n {ct,/3 = 0) is real symmetric and thus having real eigenvalues. Moreover 
its eigenvalues are simple in accordance with the following theorem (see Appendix B for 
the terminology): 

Theorem 3.2 (Trench 1993). Suppose that f is nonincreasing and 

/(0+) = M > m = /(tt-). (34) 

Then for every n the matrix A n has n distinct eigenvalues in (m,M), its even and odd 
spectra are interlaced, and its largest eigenvalue is even. 

In our case the symbol of the truncated Toeplitz matrix is a decreasing function of 9 G 
(— 7r, 0) and /(— 7T+) > /(0— ) thus the theorem is applicable. 

Proposition 3.3. The matrix A n is strictly positive definite. 

This is a consequence of the previous theorem since all the eigenvalues belong to the 
interval (/(0-), /(-7T+)) with /(0-) > 0. 

Diagonalizing A n by an orthogonal matrix and using the cyclicity property of the trace 
we obtain 

Tr (A0 = - | ^l^] = j *° j =► -«,!(«) = - f> > 0, (35) 

Icos^^Jle" |cos( z yj|e n f^i 

where A« are the eigenvalues of *A n . Also combining fl35l) and fl20l) we can express u>q in 
terms of Aj's 

/ \ (1 -oQg-i / x (1 -a)r(g-q) 1 ^ 

= T^Tiy = r(2-a)r( g + i) n ^ Ai - (36) 
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The determinant of A n is 



"*^H 2 i4?)k ) H A - >0 (37) 



since A, > 0, Vi = 1, • • • , n. 



(ii) For /3 G { — 1, 1} the operator A(a,/3 = ±1) is one-sided and the corresponding matrix is 
a sum of a lower (or upper)-diagonal matrix plus a matrix J n with elements 5^+1 = 2u>o 
(or = 2w ). For n = 2 we have two real eigenvalues while Vra > 3 we have complex 
and real eigenvalues. 

In the most general case (a G (1,2) and (3 G (—1,1)) the matrix A n (a,j3) decomposes into 
a sum of a symmetric a-dependent matrix and an antisymmetric /^-dependent matrix. The 
eigenvalues in this case are both real and complex. 



4 Asymptotic behaviour of the product of eigenvalues 
for A n 

In this section we examine the validity of the conditions under which the Szego's strong limit 
theorem holds [12] and determine the asymptotic behaviour of the product of eigenvalues. 

Definition 4.1. The set of all functions 

/(*) = z = e ie eT, with ^|/ fc |<oo (38) 

fcez fcez 

is denoted by W := W(T) and is called the Wiener algebra. 

Definition 4.2. The set of all functions ( fffff)) which belong to L 2 and satisfy 

J2(\k\ + l)\fk\ 2 <™ (39) 

fcez 

is denoted by B^ 2 and is referred as a Besov space. 

1 /2 

Proposition 4.3. The symbol f given by (27\ ) belongs to the space WnB^, it has a zero on 



T at 9 = and the winding number of f around the origin vanishes, wind(f, 0) = 0. 
Proof. □ 
/ is an element of the Wiener algebra since 

£l/»l = 2 | co f( ?) | (\M + lt(\f-n\ + \fn\) 



n=l 

oo 



l cos (-)l V J l cos (- 



2 



8 



Also / belongs to the Besov space since /ei 2 and (see (1A.4K and (1A.6j) for the proof) 

K 2 

%y (4^ + lb(|w | + |w 2 |r + S 

n=2 



£(H + l)l/n| 



< 



4 cos 2 (2p) 



4w 2 + 16(|w | + |w 2 |) 2 + 8^ (n + l)w n+l 



COS' 



\ 2 ) 



a ( -a — 1 

4 



16 

9tt 



< oo. 



(41) 



The symbol / has a zero on T at 9 = 0. Transversing T once the argument of f(e ie ) will return 
to its original zero value and the winding number vanish. 

1 /2 

Theorem 4.4 (Szego's strong limit theorem). If f G W PI _B 2 /ias no zeros on T and 
wind(f, 0) = then 



where 



ra->oo ^\J) 



G(f) = exp(log/) 

(oo 
j>(log/)*(log/)- 
fe=i 



(42) 

(43) 
(44) 



and (logf)k the Fourier coefficients o/log(/(e )). 

In our case all requirements are fullfiled except that / has zeros at 9 = 2rmr, m G Z. 
Subtracting these points from T and applying the theorem for a G (1,2) and (3 G [1, —1] the 
Fourier coefficients, using the cosine and sine series, are given by 



(log/)* = — / \og(f(e*))e-^d9 



1 

2^ 



log(n 2 (fl) + v 2 (9)) cos(k9)d9 + i / lo 



u{9) +iv{9) 
u(9)-iv(9) 



sin(k9)d9 



(45) 



The calculation of such integrals is cumbersome and we investigate only the S a (c) laws. Equa- 
tion (145!) is reduced to 



(log/)/ 



7T 



log(/(O)cos(fc0)d0 



(46) 



bearing in mind that \og(f (e 10 )) is an even function. 

Theorem 4.5. The asymptotic behaviour of eigenvalues for the problem / fl?3j) is given by (54\ ). 
Proof. □ 
The generating function 

,001 , 



log(/(e ie ))=logA^-log|cos( 



log 



9 


a 






2 cos - 


+ log ^COS 




) 











(47) 



is T = 27m! periodic (with 9 G (— 7T, 7r)) and using (|46|) we study the following two cases 
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(i) k=0. The corresponding Fourier coefficient is given by 



log 



cos 1 



, a7r 1 

I 71 

a 



P7T 


9 


/ log 


2 cos - 


to 


2 



d6 



de. 



The first integral on the right handside vanishes since 



A7T 


e 




6 




/ log 


2 cos - 


d0= / log 


2 sin - 


d0 = 


to 




Jo 


2 





The second one after the change of variable, 

a 



1- 



becomes 



log 



cos [0(1-- 



(16 



7T 1- 



(2 - a)7i 
2 

~(2-a)n 



L U 1 



log(cos u)du 



where L is Lobachevskiy's function [16]. Thus 



(log /) = log K a - log | cos 



air. 



2 71 (2 - a)7r 
(u) G Z+. The Fourier coefficients are given by 



L 7T 1- 



(log/)/ 



+ 



n 

7T 

1 







7T 

/l + / 2 



log 
log 



2 cos 







cos(k6)d6 



a 



cos (0(1--)) cos(k6)d6 



with 



and 



(2 - a)7T J 



K^f) 



a 



; 2k 



log(cos(w)) COS 



2ku 
2 - a 



du 



2kn 



7T 1- 



cos 



(l3| - lj U - COS (jK + lj 



cos u 
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The second equality in (1551) is justified by partial integration. Thus (T5TT) is written as 

r* 1 /"r(l-f) /COS (jAr - 1] It - COS (jAr + 1] It' 



We study two physically well-known cases: 



(1) The Gaussian law with a = 2 and (5 — 0. Substituting f(e ld ) = K 2 (2sin|) 2 into (146]) we 
obtain 



(log /) = log K 2 + - [ log (2 sin ( ? I I d6 = log A 2 . 



(i) k = 



7T./n 

(ii) 7^ and # 7^ 2m7r, m G Z 



(log/)* = ^ jT log ^2 sin cos(k9)d9 = ~, k e N. (56) 



Applying theorem (1421) we get 



det(A n (a,p = 0)) E » 1 c 
hm — = e z ^ fe=1 fc = e hm m (57) 



where the Euler's constant C is defined by C = linv^^ (J^fcLi ¥ — logm). Expressions 
(I57p and (|15|) differ only by the constant factor e c . This result is also justified by (39) of 
in the limit of the unit circle R = 1 and setting a = 0. 



(2) The Holdsmark law with a — | and ft = [TfJ [18]. The Fourier coefficients from ( 150]) and 
(154)) are found to be 

(bg/(a=|)) = logA^-log^ + ^G (58) 

\ m=l 

where G is Catalan's constant defined by 



(2m 

m=0 v 

and the ^ function is defined by \l/(x) 



G = V = 0.915965 • • • 

^„ 2m + 1 2 



diogr(x) 

da; 



When a takes values on the Farey series @ then the zero and nonzero Fourier coefficients are 
expressed in terms of polylogarithmic and hypergeometric functions. 



The Farey series T n of order 11 is the ascending series of irreducible fractions between and 1 whose 
denominators do not exceed n. Thus a — ^ belongs in T n if 

< p < q < n, (p,q) = l 
where (, ) denotes the highest common divisor of two integers. 
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5 Conclusions 



At quantum level using the infinitesimal generator of time translations with vanishing skew- 
ness and the definition of the Griinwald-Letnikov fractional derivative for the one-dimensional 
infinite well potential with Dirichlet boundary conditions, we established the correspondence 
between the disctretized version of the fractional Schrodinger problem and Toeplitz matrices. 
The next step was to check the conditions under which the Szego's strong limit theorem is valid 
and subsequently to determine the asymptotic behaviour of the product of eigenvalues given 
by ( )44|) and ( 154|) . An open question which deserves investigation is the general S a (j3,c,r) law. 

Finally, from the physical point of view there is a plethora of applications related to the 
fractional diffusion equation [10] but not to the Schrodinger equation. Recently, the authors 
of [19] constructed a one- dimensional lattice model with a hopping particle and numerically 
obtained the eigenvalues and eigenfunctions in a bounded domain with different boundary con- 
ditions. It would be interesting to study the asymptotic behaviour of the product of eigenvalues 
in such a model using the Griinwald-Letnikov fractional derivative. 



Appendix A 

Proof of (EU) 

The infinite series (1211) is written as 



em>" ( : ) 

n=0 v 7 



a(a — l)T(—a) 



2n, 



1" 



—a) 



T(2-a) ^ ~' v Wn l> + 1) 
= 1 F (-a;l) = 0, ae(0,2)/{l} 
and converges to zero since iF (—a; z) = (1 — z) a . 



(A.1) 



Proof of (1281) 

We prove first the identity 



E 



D T ^ + 1 „ ^ = -r(2-a) (e ie 2 Fi{l,2-a;3;e ie ) + e^F^l, 2 - a; 3; e"* 9 )) 



n=l 

ir(-a) 



T(n + 2) 4 



2 a (sin- cos - (9(2 - a) + vra) + 2(a - cos 9) 



Writing cos(9) = \(e + e ) it is enough to calculate 

» r(n + l-q) in6 ie T(2-a) (1) B (2 - a) n e md 



i r(n + 2) 



r ( 3 ) ( 3 )« 



m 



r(3) 

e~ ld Y(-a) [(I - e ie ) a - 1 + ae w ] 
r(-a) [iF (-a;e ie ) - 1 + ae 



-i6i 



iO 



where (a) n is the Pochhammer's symbol 



(a) n = a(a — 1) ■ ■ ■ (a + n — 1) 



T(a + n) 
Via) 



a) = 1, (l) n = n! 



(A.2) 



(A.3) 
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The second term of the sum (1A.2j) is derived from (IA.30 by complex conjugation. 

Proof of (Hljl 

In (HU) we used that the sum 



1 Ira + 1 — a 



n=2 n=2 

converges. It is easily checked using D'Alembert's test 



< 1, Vn > 2 (A.5) 



a n+ i _ (n + 1) (n + 1 - a) 
a n n (n + 2) 2 

and the positivity of a n 's. Also we used the identity 

g( r( n + 2) ) -("Wv ,f ' M "°' J "" ;4 * ) ' (A ' 6) 
Appendix B 

Let 

Jij := 5 ijn+1 ^j, i,j = 1,2, • • • ,n (B.l) 

be the matrix with ones along the secondary diagonal and zeroes elsewhere. For a real, sym- 
metric and Toeplitz matrix one can prove that [13] 



[J,A n }=0 (B.2) 

Using ( IB. 21) and since J 2 = I this implies that ^4 n x = Ax if and only if ^4 n (Jx) = A(Jx). If A 
has multiplicity one then from || x H2HI ^ x II 2 an d ^ = cx we conclude that c = ±1. 

Definition 5.1. An n- dimensional vector x will be called symmetric if 

Jx = x (B.3) 

or skew- symmetric if 

Jx = -x. (B.4) 

Definition 5.2. The eigenvalue A of A n is even (or odd) if it has an associated symmetric (or 
skew-symmetric) eigenvector x. 
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